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Abstract Wc investigate least-squares fitting methods for estimating the wind- 
ing rate of field lines about the axis of twisted magnetic-flux tubes. These meth- 
ods estimate the winding rate by finding the values for a set of parameters that 
correspond to the minimum of the discrepancy between vector magnetic-field 
measurements and predictions from a twisted flux-tube model. For the flux-tube 
model used in the fitting, we assume that the magnetic field is static, axisym- 
metric, and does not vary in the vertical direction. Using error-free, synthetic 
vector magnetic-field data constructed with models for twisted magnetic-flux 
tubes, we test the efficacy of fitting methods at recovering the true winding 
rate. Furthermore, we demonstrate how assumptions built into the flux-tube 
models used for the fitting influence the accuracy of the winding-rate estimates. 
We identify the radial variation of the winding rate within the flux tube as one 
assumption that can have a significant impact on the winding-rate estimates. We 
show that the errors caused by making a fixed, incorrect assumption about the 
radial variation of the winding rate can be largely avoided by fitting directly for 
the radial variation of the winding rate. Other assumptions that we investigate 
include the lack of variation of the field in the azimuthal and vertical directions in 
the magnetic-flux tube model used for the fitting, and the inclination, curvature, 
and location of the flux-tube axis. When the observed magnetic field deviates 
substantially from the flux-tube model used for the fitting, we find that the 
winding-rate estimates can be unreliable. We conclude that the magnetic-flux 
tube models used in this investigation are probably too simple to yield reliable 
estimates for the winding rate of the field lines in solar magnetic structures 
in general, unless additional information is available to justify the choice of 
flux-tube model used for the fitting. 
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1. Introduction 

Methods for reliably estimating the twist present in solar magnetic-flux concen- 
trations are required for several areas of solar-physics research. Some examples 
include: The latitudinal variation of the twist present in active regions (e.g. 
Pevtsov, Canfield, and Metcalf, 1995; Longcope, Fisher, and Pevtsov, 1998; 
Pevtsov, Canfield, and Latushko, 2001; Nandy, 2006) provides constraints for 
models of the generation and dynamics of magnetic-flux tubes in the solar 
convection zone (e.g. Rust, 1994; Linton, Longcope, and Fisher, 1996; Longcope 
and Fisher, 1996; Longcope, Fisher, and Arendt, 1996; Longcope, Fisher, and 
Pevtsov, 1998; Fan, Zweibel, and Lantz, 1998; Oilman and Charbonneau, 1999; 
Abbett, Fisher, and Fan, 2000; Choudhuri, Chatterjee, and Nandy, 2004). The 
twist present in a magnetic-flux tube can provide information about the suscep- 
tibility of the flux tube to the kink instability, where magnetic twist is rapidly 
converted to writhe (deformation of the flux-tube axis), which may play a role 
in eruptive events such as coronal mass ejections (e.g. Raadu, 1972; Hood and 
Priest, 1979; Craig et al., 1990; Vclli, Hood, and Einaudi, 1990; Linton, Longcope, 
and Fisher, 1996; Lionello et al, 1998; van dcr Linden and Hood, 1998, 1999; 
Baty et al, 1998; Baty, 1997, 2001; Fan and Gibson, 2003, 2004; Leamon et al., 
2003; Torok and Kliem, 2003; T6r6k, Kliem, and Titov, 2004; Leka, Fan, and 
Barnes, 2005; Nandy et al., 2008). Estimates for the twist present in photospheric 
magnetic-flux concentrations are useful inputs for coronal magnetic-field simula- 
tions (e.g. Mackay and van Ballegooijen, 2006, 2006; Nandy et at, 2008; Yeates, 
Mackay, and van Ballegooijen, 2008; Yeates and Mackay, 2009; Yeates et at, 
2010) and provide constraints for simulations of emerging twisted magnetic- 
flux tubes (e.g. Emonet and Moreno-Insertis, 1998; Hughes, Falle, and Joarder, 
1998; Fan and Gibson, 2003, 2004; Murray and Hood, 2008; Hood et al, 2009; 
MacTaggart and Hood, 2009; Cheung et al, 2010; Toriumi and Yokoyama, 2011). 
The magnetic twist is closely related to the twist component of the magnetic 
helicity (e.g. Berger and Field, 1984 Moffatt and Ricca, 1992). 

In a simple magnetic-flux tube, the magnetic twist can be decomposed into 
the product of two quantities: the length of the flux tube [I], and the winding rate 
[q] (i.e. the angle that a field line rotates about the flux-tube axis per unit length 
along the axis), assuming that the winding rate does not vary in the direction 
along the flux-tube axis. The focus of this article is methods for estimating the 
winding rate of the field lines in magnetic-flux tubes from vector magnetic-field 
measurements at a single observation height. Several such methods currently 
exist that we now briefly review. 

The ctbest method (e.g. Pevtsov, Canfield, and Metcalf, 1995; Longcope, Fisher ,| 
and Pevtsov, 1998; Leka, 1999; Leka and Skumanich, 1999; Pevtsov, Canfield, 
and Latushko, 2001; Leamon et al, 2003; Holder et at, 2004; Leka, Fan, and 
Barnes, 2005; Nandy, 2006) computes a single (constant) value of a, the force- 
free parameter in the equation VxB = aB, that provides the least-squares 
best-fit between a linear force-free magnetic-field model (computed using the 
observed vertical component of the magnetic field) and the observed horizon- 
tal components of the magnetic field, usually for an entire active region. This 
method was used by Pevtsov, Canfield, and Metcalf (1995), Longcope, Fisher, 
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and Pevtsov (1998), Pevtsov, Canfield, and Latushko (2001), Holder et al. (2004), 
and Nandy (2006) to show how abest in active regions depends on latitude. By 
averaging over the entire active region this technique is less susceptible to noise 
than methods that compute a by differentiating the magnetic field (discussed 
later). However, the large-scale averaging means that the magnetic twist is not 
probed on spatial scales smaller than an entire active region (e.g. Leka and 
Skumanich, 1999; Leka, 1999; Leka, Fan, and Barnes, 2005). Another limitation 
of the abest method is that the magnetic field at the solar photosphere, where the 
vector fields are typically measured, is not force-free (e.g. Metcalf et al., 1995). 

The abest method was also used by Leamon et al. (2003) to compute the 
winding rate and, subsequently, the twist present in active regions associated 
with eruptive events, in order to determine whether the kink instability was the 
mechanism responsible for the eruption. To determine the winding rate, Leamon 
et al. assumed a^est = 2g (see also Nandy, 2006); a similar assumption is invoked 
by Tiwari et al. (2009) and Tiwari, Venkatakrishnan, and Sankarasubramanian 
(2009) but for a method different from the ab os t method. Leka, Fan, and Barnes 
(2005) pointed out that this relationship between the force-free parameter [a] 
and the winding rate [q] is generally only true at the axis of an axisymmetric 
magnetic-flux tube, for a flux tube that is not thin with an axis that is straight 
and vertical (e.g. Ferriz-Mas and Schiissler, 1989). Furthermore, Leka, Fan, and 
Barnes (2005) showed that the abest method tends to underestimate the winding 
rate on the small spatial scales that may be relevant to determining whether a 
magnetic-flux system in an active region is susceptible to the kink instability 
(see also, e.g., Leka and Skumanich, 1999; Leka, 1999). 

Leka, Fan, and Barnes (2005) introduced the a pca k method for estimating the 
winding rate of the magnetic-field lines in the vicinity of the axis of a twisted 
magnetic-flux tube. The a pea k method assumes that the axis corresponds to the 
region of largest a z for a flux tube with a uniform winding rate. The relationship 
a z — 2q is then used to estimate the winding rate at the flux-tube axis. As 
shown in Table 1 of Leka, Fan, and Barnes the relative discrepancy between 
the inferred and true winding rates at the axis recovered with this method 
is roughly 10%, when tested with synthetic data from the simulation results 
of Fan and Gibson (2003, 2004); note that <7 pca k = a pca k/2 = —10 radians 
L~ x for the simulation results of Fan and Gibson (2003, 2004), where L is a 
length scale. Two aspects of this method make the results potentially sensitive 
to noise in the measurements. First, this method samples only a handful of 
points in the vicinity of the inferred location of the flux-tube axis. Second, a z 
is computed by differentiating the vector magnetic-field measurements (see also, 
e.g., Pevtsov, Canfield, and Metcalf, 1994; Leka and Skumanich, 1999; Leka, 
1999; Nandy et al, 2003; Hahn et al., 2005; Su et at, 2009; Tiwari et al, 2009; 
Tiwari, Venkatakrishnan, and Sankarasubramanian, 2009), and differentiation is 
generally a noise amplifying operation (unless techniques are used to suppress 
the influence of noise) . 

Nandy et al. (2008) introduced a least-squares fitting technique for estimating 
the winding rate in magnetic-flux tubes from vector magnetic-field data. To 
estimate the winding rate in a magnetic-flux tube, the fitting method outlined 
by Nandy et al. finds the values of a set of model parameters, related to the 
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winding rate, that correspond to the minimum of the discrepancy between the 
observations and the predictions from a prescribed flux-tube model (hereafter, 
the fitting model). This type of approach has several potential advantages over 
methods that use the force-free parameter. These include: i) fitting methods do 
not assume that the magnetic field is force-free; ii) fitting methods do not require 
derivatives of the magnetic-field measurements and are, therefore, likely to be 
less sensitive to noise than methods that do involve such derivatives; iii) fitting 
methods use many observation points made within a solar magnetic-flux tube 
and are, therefore, likely to be less sensitive to noise than methods that only 
use a few points; and iv) fitting methods may be able to probe the winding rate 
throughout the flux-tube interior, not only in the vicinity of the axis. However, 
there arc still several possible sources of error in the approach outlined by Nandy 
et al. that need to be investigated. For example, Nandy et al. assumed that the 
magnetic field in the fitting model is axisymmetric and does not vary in the 
vertical direction. Moreover, Nandy et al. assumed that the winding rate inside 
the flux tube varies with radius [r] according to q(r) = qet + c/r, where the best 
fit value of qat is their estimate for the winding rate and c is another model 
parameter; this equation implies that the winding rate is infinite at r — for 

The purpose of this investigation is to demonstrate how assumptions built into 
models for twisted magnetic-flux tubes can affect estimates for the winding rate 
inferred with least-squares fitting methods. For this purpose we rely on error- 
free synthetic vector magnetic-field measurements for which the correct winding 
rate is known. The outline of this article is as follows: The details of the fitting 
model and the assumptions built into it are described in Section 2. The details 
of the fitting method and the metrics used to test its performance are described 
in Section 3. In Section 4 we examine a version of the fitting method that makes 
a fixed assumption about the radial variation of the winding rate (in a similar 
fashion to Nandy et al., 2008). We test this fitting method with synthetic data 
constructed with models for twisted magnetic-flux tubes that have a radially- 
dependent winding rate that differs from the fitting model. We find that the 
inferred winding rate is sensitive to the assumed radial variation of the winding 
rate built into the fitting model. Subsequently, in Section 5 we present and test a 
fitting method that infers the radial variation of the winding rate. We find that 
this method can generally perform better than fitting methods that have a fixed 
radial variation of the winding rate, provided that the magnetic field in the fitting 
model and that used to construct the synthetic data is static, axisymmetric, 
and does not vary in the vertical direction. In Section 6 and the Appendix, 
we determine the sensitivity of the winding-rate estimates retrieved by fitting 
methods to the other assumptions built into the fitting model, such as the neglect 
of vertical and azimuthal variations in the magnetic field. We also apply the 
fitting method to synthetic vector magnetic-field measurements generated with 
a twisted toroidal magnetic-flux tube model, which is substantially different from 
the fitting model. In this case, we find that the winding-rate estimates can be 
inaccurate. Consequently, we do not apply the fitting method to real vector 
magnetic-field data. In Section 7 we draw conclusions. 
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2. The Fitting Model 

The fitting model is a simplified model for an isolated, static, twisted magnetic- 
flux tube. We assume that error-free, ambiguity-resolved, vector magnetic-field 
measurements (B x , B y , B z ) of a magnetic- flux tube are available at a set of 
discrete locations on an x-y-plane at constant height z, taken to approximate 
a layer of constant optical depth within a limited field of view. To describe 
the fitting model we work in a cylindrical coordinate system (r, 0,z), assuming 
that the flux-tube axis is straight, parallel to the z-direction, and located at 
r = 0, where r is the radial coordinate and 9 is the azimuthal angle measured 
counterclockwise from the positive x-axis to the position vector in the x - y-plane. 
We assume that the magnetic field in the fitting model does not vary in the 9- or 
z-directions. Under these assumptions, the radial component of the field in the 
fitting model must be zero for the magnetic field to be divergence-free and for 
the radial component of the field to be finite at the flux-tube axis. For this class 
of flux tubes, the winding rate of a magnetic- field line at a given position (i.e. 
the angle that a field line rotates about the flux-tube axis per unit length along 
the axis) is given by q m (r) = Bg(r)/(rB z (r)), where Bg = —sm9B x +cos9B y 
is the azimuthal component of the field (e.g. Priest, 1984). In these flux-tube 
models, the winding rate is independent of 9 and z, but can vary in the radial 
direction. 

To understand the radial variation of the magnetic field and winding rate 
permitted in the fitting model we refer to the results of Ferriz-Mas and Schiisslcr 
(1989). We assume that the magnetic field within the fitting model is finite and 
continuous, but may have a current sheet at the external boundary of the flux 
tube. At the axis of the fitting model we assume that the vertical component of 
the magnetic field is non-zero, and that the azimuthal component of the field is 
zero (by definition) . We assume that Bg and B z can both be expanded as a Taylor 
scries in the radial direction about the flux-tube axis. Under these assumptions, 
Ferriz-Mas and Schiissler showed that the Taylor-series expansion for B z can only 
contain even powers of r, whereas the expansion for Bg can contain only odd 
powers of r. It can then be shown that the Taylor-series expansion of the winding 
rate q m in the radial direction about the flux-tube axis can contain only even 
powers of r, and q m must be finite at the flux-tube axis. We note that it is possible 
to construct a combination of Bg and B Zl satisfying the above assumptions, 
such that q m has a singularity in the flux-tube interior (e.g. at locations where 
B z = 0). Thus, we restrict our attention to cases with Bg and B z such that 
the winding rate is finite throughout the flux-tube interior (see Figure 1(a) for a 
selection of example winding-rate profiles) . In this class of flux tubes, the current 
density has the following properties: its radial component is zero, the Taylor- 
series expansion for its azimuthal (vertical) component contains only odd (even) 
powers of r, and its azimuthal component vanishes at the flux-tube axis. The 
Lorentz force only has a radial component, its Taylor-series expansion contains 
only odd powers of r, and it vanishes at the flux-tube axis. 

Due to the assumptions outlined above, these fitting models may be too simple 
to provide a realistic representation of solar magnetic-flux tubes in general. 
Nevertheless, it is reasonable to assume that these models could be a useful 
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Figure 1. (a) Winding rate [q m = Bg/(rB z )] as a function of dimensionless distance 
from the flux-tube axis [r/R] where R is the radius of the flux tube. The purple curve 
is for a flux tube with a uniform winding rate, qm{r) = q c . The blue curve is for a flux 
tube with a winding rate that decreases monotonically to zero at the flux-tube boundary, 
Qm(f) = 9 C (1 _ { r /R) 2 ) 2 - The red curve is for a flux tube with a winding rate that also 
decreases monotonically, but does not vanish at the boundary, <? m (r) = q c exp(— (r/R) 2 ). The 
green curve is for a flux tube with a winding rate that is non-monotonic and vanishes at the 
boundary, q m {r) = 9 C (1 + 2(r/i?) 2 — 3(r/R) 4 ). In all cases, q = 1 radian Mm -1 and, thus, 
the winding rate at the flux-tube axis IS Qrn (0) = 1 radian Mm -1 , (b) Same as (a) except for 
the corresponding ratio, w m = Bg/B z , as a function of r/R for a flux tube with a radius of 
R = 2 Mm. 

starting point for developing methods to estimate the winding rate of the mag- 
netic field in solar magnetic-flux tubes; indeed, a similar model has been used 
by Nandy et al. (2008) to estimate twist in solar active regions. We use these 
fitting models because they are simple, generally current-carrying, generally have 
a non-zero Lorentz force (which is expected near the solar photosphere, e.g. 
Metcalf et al., 1995), and the winding rate can vary in the radial direction. 
Furthermore, methods based on these simple flux-tube models can be tested 
and the influence of the assumptions built into the fitting model can be easily 
demonstrated. One such assumption that we examine is the radial variation of 
the winding rate used in the fitting model. These types of magnetic flux-tube 
models have been used in a couple of different solar contexts to understand the 
influence of different functional forms for the radial dependence of the winding 
rate. For example, they have been used in theoretical studies to determine how 
the radial variation of the winding rate affects the susceptibility of these types 
of flux tubes to the kink instability (e.g. Raadu, 1972; Hood and Priest, 1979; 
Einaudi and van Hoven, 1983; Craig et al., 1990; Mikic, Schnack, and van Hoven, 
1990; Velli, Hood, and Einaudi, 1990; Lionello et al., 1998; van der Linden and 
Hood, 1998, 1999; Baty et al, 1998; Baty, 1997, 2001). They have also been 
used in numerical simulations of rising flux tubes to understand how the radial 
variation of the winding rate influences the emergence process (e.g. Emonet 
and Moreno-Insertis, 1998; Hughes, Falle, and Joarder, 1998; Murray and Hood, 
2008). 
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3. Estimating the Winding Rate in a Magnetic-Flux Tube 

To estimate the winding rate in a magnetic-flux tube we follow Nandy et al. 
(2008) by formulating a classical fitting problem that seeks to minimize the 
discrepancy between the observed ratio w = Bg/B z and that produced by the 
fitting model, 



where r, is the radial coordinate of observation point i, N is the number of mea- 
surements used for the fitting, w b s ,i is the ratio Bg / B z derived from the vector 
magnetic-field measurements at observation point i, and w m (r,i) = riq m (ri) is 
the value of the ratio predicted by the fitting model at rj. For the fitting model, 
the Taylor-series expansion of w m in the radial direction about the flux-tube axis 
can contain only odd powers of r, and w m must vanish at the flux-tube axis (for 
some examples of w m (r) see Figure 1(b)). 

The definition of x 2 given by Equation (1) could be modified in standard 
ways to include measurement uncertainties if they were available. We use this 
definition because we are assuming (for demonstration purposes) that the mea- 
surements are free of error. In practice the ratio w — Bg / B z may not be the 
best quantity for fitting because B z is in the denominator and this may amplify 
noise and/or uncertainties in the measurements. As mentioned previously, to 
avoid situations where w is undefined, we assume that B z is non-zero at all of 
the observation points for both the fitting model and the measurements. We use 
the ratio w = Bg/B z for the following reasons: it is sufficient for demonstrating 
our main findings (with error- free measurements), it was employed by Nandy 
et al. (2008), and it does not involve directly fitting for both components of 
the magnetic field, Bg and B z (i.e. a model for each of the components of the 
magnetic field is not required). 

In what follows, we will determine the sensitivity of the winding rates esti- 
mated by this method to the assumptions built into the fitting model, using 
synthetic vector magnetic-field measurements generated with models for twisted 
magnetic-flux tubes. In general, the objective of fitting a model to data in this 
way is to extract an underlying trend from the data and possibly to estimate 
some parameters related to the best-fit model. To test the efficacy of the fitting 
method we focus on two quantities: the winding rate at the flux-tube axis and 
the radial average of the winding rate, 



where [a, b] is the range over which the average is taken. We use the on-axis and 
average winding rates as metrics because they can be used to characterise the 
twist present in a magnetic-flux tube and both of these quantities may be useful 
for determining whether a flux tube is unstable to the kink instability (e.g. Mikic, 
Schnack, and van Hoven, 1990; Velli, Hood, and Einaudi, 1990; Robertson, Hood, 




(1) 




(2) 
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and Lothian, 1992; Baty, 2001). We note that if sufficient measurements of w b s ,i 
are available then these quantities may be approximated directly from the data 
(depending on the distribution of the observation points r,) and, therefore, the 
fitting procedure may not be strictly necessary for error-free data. However, we 
emphasise that the goal here is to examine how these quantities, as estimated 
from the best-fit model, are affected by the assumptions built into the fitting 
model. 



4. Fitting Models with a Fixed Radial Dependence for the Winding 
Rate 

To demonstrate how assumptions about the radial variation of the winding rate 
for the fitting model can affect estimates for the on-axis and average winding 
rates, we consider a class of fitting models with q m (r) — q c f(r), where q a is 
a parameter with dimensions of a winding rate (i.e. radians per unit length) 
and f(r) is a dimensionless function of r that describes the radial variation of 
the winding rate (and satisfies the assumptions described in Section 2, i.e. the 
Taylor-series expansion of / about r = contains only even powers of r). In 
this section, we assume that the functional form of / is fixed during the fitting 
procedure and that / is independent of q c (i.e. df /dq c = 0). Consequently, the 
minimization problem described by Equation (1) has only one free parameter 
[q c ] that controls the amplitude of the winding rate for the fitting model. 

By fixing the functional form of / we are making an assumption about the 
radial variation of the winding rate. This may be reasonable if observational or 
theoretical information is available to support this assumption, and it may be 
a useful starting point for developing fitting methods to estimate twist (e.g. 
Nandy et al., 2008). To understand the implications of this assumption we 
construct error-free synthetic data using flux-tube models that satisfy the same 
assumptions as the fitting model (i.e. the magnetic field and winding rate do 
not vary in time or in the azimuthal or vertical directions, see Section 2 for 
details). To this end, in the minimization problem described by Equation (1) 
we set w bs,i — q c obs r ig( r i), where q c obs is a parameter with dimensions of a 
winding rate and g(r) is a dimensionless function of r that satisfies the same 
assumptions as f(r). Thus, for this exercise the only differences between the 
model used to construct the synthetic data and the fitting model are the values 
of q a and q a oba , and the radial dependence of the winding rates described by / 
and g. To generate synthetic data we sample w b s .i along a line at fixed 0, which 
is appropriate because none of the quantities of interest vary in the azimuthal 
direction. 

For this combination of synthetic data and fitting model, the minimization 
problem can be solved exactly. It can be shown that the best-fit value for q c , 
corresponding to the minimum of % 2 defined by Equation (1), is q* = c w q c obe , 
where the asterisk denotes the best-fit value and 



(3) 
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Table 1. Results for tests of a fitting model with a uniform winding 
rate [q m (r ) = q c ] applied to several different synthetic data sets with 

"obs,. = q c , ohe ng{n). 



Synthetic data generated with g(r): 


9o/<?0,obs 


^av/^av^bs 


x 2 /(« 2 < obs ) 


g(r) = (1 - (r/Rff 


0.227 


0.425 


0.0196 


g(r) = 1 + 2(r/i?) 2 - 3(r/R) 4 


0.907 


0.85f 


0.0557 


g(r) = c x p(~(r/Rf) 


0.567 


0.759 


0.00806 



Consequently, the winding rate for the best-fit model is q* (r) = c w q c ohs f(r). For 
synthetic data with g(0) ^ 0, the inferred value for the winding rate at the flux- 
tube axis is directly proportional to the true value, i.e. <|q = qo, bsC w f{0)/g(0), 
where the in the subscript denotes the value at the axis. The inferred and true 
average winding rates are also directly proportional, i.e. 

* _ Jaf( r ) dr 

<7av — <7av,obsCto— j — — - , (4) 

J a 9{r)dr 

for synthetic data with f g(r)dr ^ 0. In each case, the constant of proportion- 
ality relating the inferred and true winding rates is determined by the mismatch 
between the radial variation of the winding rate assumed for the fitting model / 
and that used to construct the synthetic data g, at the set of observation points 
Ti and at the axis r = or over the averaging interval [a, b]. However, it should be 
noted that, for cases with / ^ g at some of the relevant locations, the constant 
of proportionality can theoretically be unity (meaning that the inferred winding 
rate matches the true value). 

To demonstrate the implications of these results, we test a fitting model with 
a uniform winding rate, i.e. f(r) = 1. A similar fitting model was used by 
Nandy et al. (2008). The difference between this approach and that taken by 
Nandy et al. is that their fitting model has two free parameters, and they did not 
require w = at r = 0. We construct three different synthetic data sets using 
different choices for the radial variation of the winding rate (see, e.g., Figure I): 
g(r) = (1 - (r/R) 2 ) 2 , g(r) = 1 + 2(r/R) 2 - 3(r/R)\ and g(r) = cxp(-(r/i?) 2 ), 
where R is the radius of the flux tube. These synthetic winding-rate profiles are 
chosen for the following reasons: they satisfy the requirement that the Taylor- 
series expansion of winding rate about r — contains only even powers of r, 
they are single signed in the interval < r < R, and profiles like these have been 
used previously by Hood and Priest (1979), Mikic, Schnack, and van Hoven 
(1990), Lionclio et al. (1998), Baty et al. (1998), and Baty (2001). Synthetic 
measurements are generated at two hundred equally spaced observation points 
in the interval < < R, all of which are used for the fitting, and the average 
winding rate is calculated for the interval [0,R]. 

The results for these tests (see Table 1) show that the fitting model with a 
uniform winding rate does not retrieve the true value for the winding rate at 
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the flux-tube axis or the average winding rate for any of the synthetic data sets; 
that is, the ratios <7o/<Zo,obs and ql v /q av ,obs both deviate substantially from unity 
(in some cases by several tens of percent). In other experiments with different 
combinations of / and g we find similar outcomes (results not shown). For the 
best-fit model, x 2 is directly proportional to <z 2 oba , i.e. 

^^yBWtn)-^)) 2 , (5) 

and the results in Table 1 show that smaller values of x 2 do not necessarily 
correspond to more accurate winding-rate estimates. 

For the tests in Table 1 the fitting model underestimates the magnitudes of 
both the on-axis and average winding rates. In other experiments we find cases 
where the fitting model with f(r) — I overestimates the magnitude of the wind- 
ing rates. Whether a particular fitting model / underestimates or overestimates 
the magnitude of the winding rate depends on the functional forms of / and g, 
and on the locations of the observation points j-j. 

The fitting model tested in Table 1 retrieves the correct sign for both the on- 
axis and average winding rates. This is expected for this combination of synthetic 
data and fitting model, and can be understood by inspecting the solutions for the 
best-fit winding rates presented above. At the axis, the inferred and true winding 
rates will have the same sign if both / and g satisfy the following conditions: 
they are single signed for all of the observation points rj and have the same sign 
at the axis that they have at the observation points r^. Likewise, the inferred 
and true average winding rates will have the same sign if both / and g satisfy 
the following conditions: they arc single signed for all of the observation points 
Ti and in the averaging interval [a, b] , and they have the same sign in the interval 
[a, b] that they have at the observation points rj. 

5. Fitting Methods that Infer the Radial Dependence of the 
Winding Rate 

In the previous section, we demonstrated how an error can be introduced into 
the estimates for the winding rates when the radial variation of the winding 
rate assumed for the fitting model is fixed and does not match that in the flux- 
tube model used to generate the synthetic measurements. In this section we 
show how this problem can be addressed by inferring the radial variation of the 
winding rate during the fitting procedure (see, e.g., Chapter 15.4 of Press et at, 
1992). To this end, we assume that the winding rate in the fitting model can be 
approximated as 

U0 = E3A(r), (6) 

where hj(r) are dimensionless basis functions, Qj are amplitudes with dimen- 
sions of a winding rate, and rib is the number of basis functions used in the 
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approximation. The basis functions [hj] must be chosen such that the radial 
variation of the winding rate satisfies the conditions discussed in Section 2 (i.e. 
the Taylor-series expansion for the winding rate [q m ] in the radial direction about 
the flux-tube axis contains only even powers of r, and q m is finite throughout 
the flux-tube interior). 

For this approach, the aim is to find the values for the set of amplitudes 
[Qk] for a given set of basis functions [hk] that correspond to the minimum of 
X 2 (as defined by Equation (1)). The resulting system of linear equations (for 
dx 2 /dQk = 0, k = 1, . . . , n b ) is: 

m 

^ AkjQ = Ck , (7) 

where 

N N 

A kj = y] r i h k{n)hj(ri) , and C k = J~] rjhk(ri)w ohSt i , (8) 
i=i »=i 

which can be solved using standard methods of linear algebra (provided that 
the matrix Akj is not singular or ill-conditioned, see, e.g., Chapter 2 of Press 
et at, 1992); in the examples that follow we use singular value decomposition 
and monitor the condition number of the matrix Akj to ensure that it is not 
ill-conditioned. 

To test this approach we follow the same procedure that we used in Section 4. 
That is, we construct error-free synthetic data using a flux-tube model that 
satisfies the same assumptions as the fitting model (i.e. the magnetic field and 
winding rate do not vary in time or in the azimuthal or vertical directions) 
by setting w b s ,i = Q c ,obe r i9( r i)- For this type of synthetic data it is evident 
from Equations (7) - (8) that the best-fit amplitudes Qj are directly proportional 
to g Ci0bs . Therefore, the best-fit on-axis and average winding rates are directly 
proportional to the respective true values (for synthetic data with g(0) ^ 
and g(r)dr ^ 0), and \ 2 f° r the best-fit model is directly proportional to 
^cobs- For demonstrating this approach we use even-order polynomials for the 
basis functions [hj(r) — (r/i?) 2 ^ -1 )] which are a convenient choice because they 
satisfy the required properties for the Taylor-series expansion of the winding 
rate about the flux-tube axis. We use two different values for rib [nb — 2 and 
rib = 3]. We note that polynomials may not be the best choice for the basis 
functions in practice, but determining the optimal type of basis functions is 
beyond the scope of this investigation. Synthetic measurements are generated 
at two hundred equally-spaced observation points in the interval < r, < R 
using g(r) = exp(— (r/R) 2 ). We use this synthetic winding-rate profile because 
its Taylor-series expansion is infinite (i.e. it cannot be represented exactly by 
a low-order polynomial approximation). The average winding rate is calculated 
for the interval [0,R]. 

We compare the results of this approach with those from several fitting mod- 
els that assume a fixed radial dependence for the winding rate (as described 
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Table 2. Results for tests of several fitting models applied to synthetic data 
constructed with u) obs ,i = <J C >obe Ti exp(-(r;/.R) 2 ) 

Fitting model g ( */<?0,obs 9av/9av,ob B X 2 /(-R 2 ^ obs ) 



9™ 


(r) 


= ?c 


0.567 


0.759 


0.00806 


1-m 


(r) 


= (?c (l-(r/.R) 2 ) 2 


1.507 


1.076 


0.0320 


1 m 


(r) 


= <? c (1 + 2(r/R) 2 -3(r/R) 4 ) 


0.571 


0.815 


0.00764 


Vm 


(r) 


= Qi+Q2(r/R) 2 


0.920 


0.970 


0.000132 


<!,„ 


(r) 


= Qi+Q2(r/R) 2 + Q 3 (r/R) 4 


0.991 


0.998 


9.46x10" 



in Section 4). For this synthetic data set, the results for these tests (see Ta- 
ble 2) show that both the on-axis and average winding rates retrieved by the 
best-fit polynomial basis-function approximations are more accurate than those 
retrieved with the fitting models with a fixed radial variation for the winding 
rate. The corresponding x 2 values are also much smaller for the best-fit poly- 
nomial basis-function approximations. For low-order polynomial basis-function 
approximations, we find that the winding-rate estimates and corresponding % 2 
values improve as n& increases. In other experiments with the same set of basis 
functions, but different types of synthetic data, we find a similar pattern (see 
Appendix A). For fitting methods with a fixed assumption for the radial variation 
of the winding rate (as described in Section 4) , the results in Table 2 show that 
different assumptions for the radial variation of the winding rate can produce 
different estimates for the winding rates (given the same synthetic data). 



6. Testing Other Fitting-Model Assumptions 

In the previous section, we demonstrated that fitting methods that infer the 
radial variation of the winding rate can retrieve more accurate estimates for 
the on-axis and radially-averaged winding rates than methods that make a 
fixed assumption about the radial variation of the winding rate (if the assumed 
radial variation for the fitting model does not match that used to construct the 
synthetic data). Those tests used synthetic data constructed with models that 
satisfy the same assumptions as the fitting model, except for the assumption 
regarding the radial dependence of the winding rate, which was considered an 
unknown. As discussed above, the fitting models used here are very simple 
and may not provide a realistic representation for twisted solar magnetic-flux 
tubes in general. The purpose of this section is to determine whether the fitting 
methods developed up to this point, which are based on simplified models for 
twisted magnetic-flux tubes, can be expected to retrieve accurate winding-rate 
estimates for synthetic data constructed with models that do not satisfy all of 
the assumptions built into the fitting model. To re-iterate, the assumptions built 
into the fitting model include: the flux tube is static and isolated, the flux-tube 
axis is straight and parallel to the ^-direction, the location of the flux-tube axis 
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is known, the fitting model is axisymmetric (i.e. the magnetic field and winding 
rate in the fitting model do not vary in the ^-direction) , and the magnetic field 
and the winding rate do not vary in the ^-direction. In the Appendix we conduct 
a series of experiments that document the effects of some of these assumptions 
on the winding-rate estimates. For the sake of brevity, here we briefly summarise 
the main findings from each experiment, before applying the fitting method to 
synthetic data generated with a twisted toroidal magnetic-flux-tube model. 

Regarding variations of the magnetic field in the vertical direction for the 
model used to construct the synthetic data (see Appendix A for details and 
examples), if the observed flux tube is axisymmetric, with an axis that is straight 
and vertical, we find that reasonable winding-rate estimates can be retrieved for 
the plane where the magnetic field is measured by the fitting methods that 
infer the radial variation of the winding rate. The winding rate at the axis of 
an axisymmetric magnetic-flux tube, with an axis that is straight and vertical, 
does not vary along the length of the axis (e.g. Ferriz-Mas and Schiisslcr, 1989). 
Therefore, an estimate for the winding rate at the flux-tube axis retrieved using 
measurements at one height can be applied along the entire axis of the flux 
tube, provided that the observed flux tube is indeed axisymmetric, with an axis 
that is straight and vertical. On the other hand, away from the flux-tube axis, if 
the magnetic field varies with height we find that the winding rate determined 
for a given radial location (or field line) at one height may not correspond to 
the winding rate at the same radial location (or field line) at a different height. 
Thus, the average winding rate found by fitting methods using measurements 
at one height cannot generally be applied to other heights, without additional 
information. 

Regarding variations of the magnetic field in the azimuthal direction for the 
model used to construct synthetic data (see Appendix B), it is clear that fitting 
methods that use axisymmetric fitting models cannot retrieve the azimuthal 
dependence of the winding rate. With that in mind, we find that fitting methods 
that use axisymmetric fitting models can retrieve reasonable estimates for the 
azimuthally-averaged, radially-averaged winding rate and the on-axis winding 
rate for synthetic measurements constructed with a magnetic field that varies 
in both the azimuthal and radial directions (but not in the vertical direction). 
For example, we find that the relative discrepancy between the inferred and true 
winding rates is of order 10% when the winding rate varies in the azimuthal 
direction by up to a factor of three in some parts of the observed flux tube. 

In Appendix C we find that an error in the location of the flux-tube axis 
generally produces an error in the winding rates retrieved by fitting methods, 
in the absence of other sources of error. The relative error in the winding-rate 
estimates is of order 10 % when the error in the location of the flux-tube axis 
is of order 10 % of the flux-tube diameter, for the cases we have examined. We 
also find that the error in the winding-rate estimates increases as the error in 
the location of the flux-tube axis increases. 

6.1. Testing the Fitting Methods with a Twisted Toroidal Magnetic-Flux Loop 

To test the efficacy of the fitting methods applied to solar-like magnetic fields 
we employ the twisted, toroidal, magnetic-flux-loop model used to drive the 
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simulations of Fan and Gibson (2003, 2004) and Cheung et al. (2010). This 
model is useful for constructing synthetic data because it qualitatively resem- 
bles some twisted solar magnetic structures (depending on the choice of model 
parameters, see, e.g., Leka, Fan, and Barnes, 2005), the winding rate of the field 
lines about the flux-loop axis is well defined, and the model violates several of 
the assumptions built into the fitting model (this partly depends on the choice of 
model parameters). Moreover, the simulation results of Fan and Gibson (2003, 
2004) have been used for "blind" testing methods for analysing vector-field mea- 
surements, such as for twist estimation (e.g. Leka, Fan, and Barnes, 2005) and 
azimuthal ambiguity resolution (e.g. Metcalf et al., 2006; Li, Amari, and Fan, 
2007; Crouch, Barnes, and Leka, 2009). 

To describe this magnetic field we use a spherical coordinate system (p, r, (f>), 
where <j> is the angle measured counterclockwise from the positive z-axis to the 
position vector in the x-z-plane, r is the polar angle measured from the positive 
y-axis to the position vector, and p is the distance from the position (x,y,z) 
to the origin (i.e. p 2 = x 2 + y 2 + z 2 ). The flux-loop axis is circular, lies in the 
x-z-plane (see Figure 1 of Fan and Gibson 2003, 2004), and is located at r = 7r/2 
and p = R t - The magnetic field is given by 



B = Vx 



A(p,r) 
psinr 



B 4> (P> T ) e 4> ■ 



(9) 



where 



A(p,r) = -q t a 2 B t exp 



(P,r) 



B<t> (p, t) 



<hBt 
psinr 



- exp 



(P,r) 



( 10 ) 

B t is a magnetic-field strength, vo 2 (p, r) = p 2 + R 2 — 2pR t sin r is the squared 
distance from the position (x, y, z) to the flux-loop axis at constant <j>, a t is a 
length scale, and qt is a dimensionless parameter that is related to the rate at 
which the field lines wind around the flux-loop axis. As in Fan and Gibson, we 
truncate the magnetic field to zero for w > 3a t . 

To understand some features of this magnetic field, it is helpful to con- 
sider a coordinate system {w,i/},<j>) that results from two cylindrical coordinate 
transformations: 



2 , 2 
X + Z , 



a sm< 



CTCOS( 



(11) 



and 



w 2 = s 2 + y 2 , w cos tp — s , msmip = y, (12) 

where w and <p are as defined above, s — Rt — <J, and ip is the angle about the 
flux-loop axis measured counterclockwise from the positive s-axis to the position 
vector in the s-y-plane at constant cf). In this coordinate system it can be shown 
that _B ro = and B^ = (q t / 'a^zuB^ (e.g. Hood et at, 2009, MacTaggart and 
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Hood, 2009). Thus, the winding rate of the field lines about the flux-loop axis, 
per unit length along the axis, is B^,/(zuB ( f > ) — qt/a t , which is constant, yet 
and Bq vary in both the w- and ^-directions but not in the ^-direction. 

To construct synthetic data we use the same parameter values used by Fan and 
Gibson (2003, 2004). We sample the magnetic field at a set of discrete locations 
on an z-y-plane at constant height zq, with a grid spacing of 0.00625 L in both 
the x- and y-directions (where L is a length scale, corresponding to the size of the 
domain in the y-direction). For the number of observation points in the x- and y- 
directions we set n x = 240 and n y = 160. We set R t = 0.375 L and a t = 0.1 L. At 
the flux-loop axis the only non-zero component of the magnetic field is parallel to 
the axis (i.e. B^, = at w = 0), and it has a magnitude of atB t / Rt- We choose 
the value of B t such that the magnetic-field strength at the axis is 2.4 kG, and 
the positive-polarity footpoint is located in the x < half-space. A typical value 
for the winding rate is chosen by making the following considerations: The twist 
T for field lines in the vicinity of the flux-loop axis, over the length of the semi- 
circle coinciding with the axis, is given by T = ql = q t l/a t , where I — irR t . Fan 
and Gibson use qt = — 1 and, thus, the twist over the semi-circular axis in that 
case is —1.875 x 2tt radians. The critical twist typically quoted for the onset 
of the kink instability is \ql\ rts 2ir radians (e.g. Raadu, 1972; Hood and Priest, 
1979; Einaudi and van Hoven, 1983; Mikic, Schnack, and van Hoven, 1990; Velli, 
Hood, and Einaudi, 1990; Lionello et ai, 1998; van der Linden and Hood, 1998, 
1999; Baty et ai, 1998; Baty, 2001; Fan and Gibson, 2003, 2004; T6r6k, Klicm, 
and Titov, 2004). To test the performance of the fitting methods for a range of 
relevant winding rates, we construct 20 synthetic vector magnetograms with 20 
equally-spaced values for q t in the range — 1 < q t < 1; the case where the winding 
rate is zero is not considered. For each synthetic magnctogram we analyse the 
positive-polarity footpoint of the flux loop as observed on the plane z = z . 
For the fitting models we set R = 3a t (i.e. the radius of the flux tube). At each 
observation point with B z > 100 G and r < R wc compute the ratio w = Bg/ B z , 
using the known location of the flux-loop axis on the z = z plane to determine 
the radial coordinate r, the azimuthal angle 9, and Bg. We calculate the average 
winding rate over the interval used for the fitting (i.e. < r < max(rj)). We 
note that by using the known values for the radius of the flux loop and the 
location of the flux-loop axis we are intentionally giving the fitting methods the 
best possible chance for success, although in practice these parameters are not 
known. 

Figure 2 shows results when the magnetic field is sampled on the plane z = 0. 
This plane is special in the sense that the cylindrical directions [r, 9, and z\ are 
locally equivalent to the directions w, ip, and (j), respectively (i.e. the axis of 
the flux loop is parallel to the z-direction for z =0). Thus, on this plane the 
radial component of the magnetic field [B r ] is zero. On the other hand, Bg and 
B z are generally non-zero and vary in both the 9- and r-directions; that is, for a 
given value of r the inferred values of Bg and B z are not single-valued (evident 
as the scatter in Figure 2(b)). Hence, the assumption that the magnetic field 
is axisymmetric is violated. However, the inferred ratio w is single-valued for a 
given value of r on the zq = plane. The fitting model that assumes a uniform 
winding rate profile can therefore retrieve the correct values for both the on-axis 
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Figure 2. (a) Vector magnetic field at zq = for the toroidal flux loop with (qt/at) = —10 ra- 
dians L . Positive (negative) vertical magnetic field [B z ] is indicated by red (blue) contours 
at 100, 500, 1000, 2000 G. Horizontal magnetic field is plotted at every eighth pixel, with 
magnitude proportional to arrow length (the maximum horizontal field strength shown is 
approximately 1286 G); arrows with a horizontal magnetic-field strength less than 25 G 
are not shown, (b) Inferred cylindrical components of the magnetic field [B r , Bg, B z ] as a 
function of radial distance from the flux-loop axis [r], for the positive-polarity footpoint of 
the magnetic field shown in (a). The red, blue, and green points correspond to B z , Bg, and 
B r , respectively (note: B r = for zq = 0). (c) The ratio, w = Bg/B z , as a function of 
radial distance from the flux- loop axis [r]. The black points are the measurements derived 
from the synthetic data in (a) and (b). The purple curve is the best-fit case for a fitting 
model with q m (r) = q c , the blue curve is for q m (r) = q c (l — (r/R) 2 ) 2 , the green curve is for 
q m (r) = q c (l + 2(r/R) 2 - 3(r/R) 4 ), the gray curve is for q m (r) = Qi + Q2(r/R) 2 , and the 
orange curve is for q m (r) = Qi +Q2 (r/R) 2 + Qs(r / R) 4 ; note that the purple, gray, and orange 
curves are effectively indistinguishable here but can be discerned in subsequent figures, (d) The 
ratio of the inferred and true winding rates at the flux-tube axis [gj/ffo.obs] as a function of 
the true winding rate at the axis [f?o,oba] for the fitting models in (c). In these tests we use 
the toroidal flux-loop model of Fan and Gibson (2003, 2004) to generate twenty synthetic data 
sets with a range of winding rates qt/o-t (see text). Each point represents the result from a 
single fitting experiment. Note that the curves joining the points are only included as a guide 
and that the purple, gray, and orange curves are effectively indistinguishable here but can be 
discerned in subsequent figures. The colors are the same as in (c). (e) Same as (d) except 
for the corresponding ratio of the inferred and true average winding rates [gJv/Qav obs] as a 
function of the true average winding rate [g av ,obs]i with the average taken over the range used 
for the fitting, (f) The corresponding value of yj? ( see Equation (1)) as a function of qo,ohs- 

and average winding rates to a high degree of accuracy, despite some of the 
fitting-model assumptions being violated. The same is also true for the best-fit 
polynomial basis- function approximations. On the other hand, the fitting models 
that assume a fixed, non-uniform winding rate do not retrieve the correct wind- 
ing rates. These results indicate that the winding rates retrieved by the fitting 
methods are not sensitive to variations of the field in the 0-direction, provided 
that the winding rate itself does not vary in the ^-direction; we acknowledge that 
this situation may be unlikely in practice. These results also indicate that the 
winding-rate estimates are not strongly sensitive to the curvature of the flux- 
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loop axis on the plane zq — 0; in experiments with smaller values of Rt we find 
results similar to those in Figure 2 (results not shown) . 
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Figure 3. Same as Figure 2, except that the magnetic field is sampled on the plane 
zq = 0.14375 L. (a) Vector magnetic field at zq = 0.14375 L for the toroidal flux loop with 
(qt/at) = —10 radians L~ 1 . The maximum horizontal field strength shown is approximately 
1730 G. (b) Inferred cylindrical components of the magnetic field [B r , Bg, B z ] as a function of 
radial distance from the flux-loop axis [r], for the positive-polarity footpoint of the magnetic 
field shown in (a), (c) The ratio [w = Bg/B z ] as a function of radial distance from the flux-loop 
axis [r]. (d) The ratio of the inferred and true winding rates at the flux-tube axis [go/9o,obs] 
as a function of the true winding rate at the axis [90, obs] f° r the fitting models in (c). (e) 
Same as (d) except for the corresponding ratio of the inferred and true average winding rates 
[llv /lav obs] as a function of the true average winding rate [g aVj obs]j with the average taken 
over the range used for the fitting, (f) The corresponding value of yj? ( see Equation (1)) as 
a function of go obs- 

Figure 3 shows results for the case where the magnetic field is sampled on the 
plane z Q — 0.14375 L. Compared to the case shown in Figure 2, this value of zq 
may be taken to represent the flux loop at an earlier stage of emergence (e.g. 
Fan and Gibson 2003, 2004). On the plane z Q = 0.14375 L several of the fitting- 
model assumptions are violated, these include: the axis of the flux loop is curved 
and not parallel to the z-direction, the magnetic field varies in the z-direction 
(although the variation is quite different from that discussed in Appendix A), 
the radial component of the magnetic field is nonzero and not single- valued (see 
Figure 3(b)), and the ratio w = Bg/B z is not single- valued for a fixed value of r 
(Figure 3(c)). The latter two features are observable and indicate that the fitting 
model is not strictly appropriate for this data set. However, judging from the 
vertical component of the field B z in Figure 3(a), the contour for B z = 100 G is 
fairly circular and it may seem reasonable to apply the fitting methods without 
modification. Taking this approach, we find that none of the fitting methods 
retrieve the correct values for the on-axis or average winding rates (Figure 3(d) 
and (e)), although in some cases the estimates are within approximately 10% 
of the true values, such as for the fitting model with a uniform winding rate. 
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We also find cases where the estimates retrieved by two different fitting models 
approximately agree, yet both estimates do not agree with the true winding rate. 
The quality of the fits is universally poor (Figure 3(f)); that is, according to x 2 
the various best-fit models are effectively indistinguishable, yet they produce a 
variety of winding-rate estimates. For these test cases we find that the inferred 
winding rates are not directly proportional to true values. We also find that 
slightly better estimates for both the on-axis and average winding rates (with 
the average taken over the fitting interval) and smaller \ 2 values can be retrieved 
by restricting the fitting to observation points closer to the axis, where the scatter 
in w bs is less severe, but this is not true for all values of z . 

In other experiments we find that the performance of the fitting methods gets 
progressively worse as z increases. Figure 4 shows the case where the magnetic 
field is sampled on the plane z a = 0.306250 L. In this case, the footpoints are not 
well separated, the magnetic field is obviously not axisymmetric, the flux-loop 
axis is highly inclined with respect to the z-axis, and the winding rate clearly has 
a large apparent azimuthal dependence. None of the fitting models consistently 
retrieve an accurate estimate for the winding rates. For larger values of the true 
winding rates, the accuracy of the winding-rate estimates retrieved by the best- 
fit polynomial basis-function approximations gets worse as n& increases, contrary 
to expectations (when the magnetic field does satisfy the assumptions built into 
the fitting models). These results indicate that the fitting methods that infer the 
radial variation of the winding rate do not generally retrieve reliable winding- 
rate estimates when the observed magnetic-flux tube deviates substantially from 
the fitting model. 

7. Conclusions 

We demonstrate how assumptions built into models for twisted magnetic-flux 
tubes used in least-squares fitting methods can influence estimates retrieved for 
the winding rate of the magnetic-field lines about the flux-tube axis. The fitting 
methods that we examine estimate the winding rate in a magnetic-flux tube by 
finding the value of a parameter (or set of parameters), related to the winding 
rate, that corresponds to the minimum of the discrepancy between observations 
and predictions from a flux-tube model, in a similar fashion to Nandy et at 
(2008). For the flux-tube model used in the fitting we assume that the magnetic 
field is static, axisymmetric, and does not vary in the vertical direction. Us- 
ing error-free, synthetic vector magnetic-field data constructed with models for 
twisted magnetic-flux tubes we test the performance of these fitting methods at 
recovering the true on-axis and average winding rates. We find that the accuracy 
of the winding-rate estimates retrieved with fitting methods is sensitive to the 
details of the assumptions built into the flux-tube model used for the fitting. 

We identify the radial variation of the winding rate as one assumption that 
can have a significant impact on the winding-rate estimates. For fitting methods 
that make a fixed assumption about the radial variation of the winding rate, we 
find that a significant error can be introduced into the winding-rate estimates 
if an incorrect assumption is made about the radial dependence of the winding 
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Figure 4. Same as Figure 2, except that the magnetic field is sampled on the plane 
2o = 0.306250 L. (a) Vector magnetic field at zo = 0.306250 L for the toroidal flux loop with 
(qt/ctt) = —10 radians L _1 . The maximum horizontal field strength shown is approximately 
2312 G. (b) Inferred cylindrical components of the magnetic field [B r , Bg, B z ] as a function of 
radial distance from the flux-loop axis [r], for the positive-polarity footpoint of the magnetic 
field shown in (a), (c) The ratio [w = Bg/ B z ] as a function of radial distance from the flux-loop 
axis [r]. (d) The ratio of the inferred and true winding rates at the flux-tube axis [go/<?0,obs] 
as a function of the true winding rate at the axis [?o,obs] f° r the fitting models in (c). (e) 
Same as (d) except for the corresponding ratio of the inferred and true average winding rates 
[llv/la.v oba] as a function of the true average winding rate [q av ,obs]j with the average taken 
over the range used for the fitting, (f) The corresponding value of yx^ (see Equation (1)) as 
a function of go, oba- 



rate. We show mathematically that the magnitude of the discrepancy between 
the inferred and the true winding rates is related to the mismatch between the 
radial variation of the winding rate assumed for the flux-tube model used in the 
fitting and that used to construct the synthetic data. We subsequently show that 
fitting methods that make different assumptions about the radial variation of the 
winding rate can produce different estimates for the winding rates when applied 
to the same synthetic data. We find that best-fit models with smaller values of 
the misfit function [x 2 ] do not necessarily retrieve more accurate winding-rate 
estimates. We also find that a model that assumes a uniform winding rate within 
the flux tube does not necessarily retrieve an accurate estimate for the radial 
average of the winding rate, as one may expect. We show that the errors caused 
by making an incorrect assumption about the radial dependence of the winding 
rate can be largely avoided by inferring the radial variation of the winding rate, 
using a basis-function approximation (e.g. Press et al, 1992). 

We apply the fitting methods to synthetic data generated with a twisted, 
toroidal, magnetic- flux- loop model (e.g. Fan and Gibson, 2003, 2004) which has 
some features that are more solar-like than the other synthetic data sets used 
in this investigation. We find that the correct winding rates can be recovered by 
fitting methods if the magnetic-flux loop is sampled on a special plane where the 
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flux-loop axis is parallel to the vertical direction, and the important assumptions 
built into the model used for the fitting are satisfied. If the magnetic-flux loop is 
sampled away from this special plane the flux loop can differ substantially from 
the model used for the fitting, in ways that are expected for solar observations of 
twisted magnetic-flux tubes. In these cases we find that fitting methods generally 
fail to recover the correct winding rates, with discrepancies between the inferred 
and true winding rates of up to several tens of percent. 

We therefore conclude that the models used for the fitting, based on static, 
axisymmetric, vertically translation-invariant flux tubes, are too simple to yield 
accurate estimates for the winding rate of the magnetic field in solar magnetic 
structures in general. However, if the observed magnetic field is not substantially 
different from the simple model used for the fitting (or information is available 
to justify the choice of fitting model) these methods may be useful for providing 
a rough estimate of the winding rate for the plane where the magnetic field 
is measured, provided that fitting methods that infer the radial variation of 
the winding rate are used. Beyond rough estimates, however, a more sophisti- 
cated model for twisted magnetic-flux tubes is required for least-squares fitting 
methods to retrieve accurate winding-rate estimates in general. 
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Appendix 

A. Testing the Fitting Methods with Twisted Magnetic Fields that 
Vary in the Vertical Direction 

To show how the fitting methods perform when the observed flux tube has a 
magnetic field that varies in the z-direction, we consider an axisymmetric, linear 
force-free, magnetic field (e.g. Schatzman, 1965; Browning and Priest, 1983; 
Petrie, 2008), governed by the magnetic- flux function 



where J is the Bessel function of the first kind, m is a parameter, k 2 = m 2 + 
a 2 , and a is the (constant) force-free parameter. We assume that m is real 
and positive so the magnetic-flux function described by Equation (13) decays 
exponentially with increasing height z. The magnetic-field lines lie in the surfaces 
of constant We assume that the external boundary of the flux tube is made 
up of the field lines on the flux surface 'J = ^>o, and R is the radius of the flux 
tube at z = [^{R, 0) = ^o]- The magnetic-field components are: 



*(r,z) = $> rJi{kr)exp(-mz)/(RJ 1 (kR)) , 



(13) 



B r (r,z) = - 



Id* 

r dz 



= if> mJi(kr) exp(-mz)/(RJi(kR)) , 



(14) 
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a 

B g (r,z) = = ^ aJ 1 (kr)exp(-mz)/(RJ 1 (kR)) , (15) 
r 

B z (r,z) = ^ =* kJ {kr)exp(-™z)/(RJ 1 {kR)). (16) 
r or 

This magnetic field is useful for constructing synthetic data to test the fitting 
methods because it has some features that are consistent with the fitting model 
(i.e. the magnetic field is axisymmetric, and the flux-tube axis is straight and 
vertical). This allows us to demonstrate the influence of the various features 
of this magnetic field that are not included in the fitting model, which are: 
the magnetic field varies in both the vertical and radial directions, the radial 
component of the magnetic field is non-zero, and the radius of the flux tube 
varies in the vertical direction (i.e. the radial location of the flux surface with 
* = $0 varies in the z-direction) . 

For a/0 this magnetic field is twisted (i.e. Bg ^ 0, see Equation (15)). At a 
given position, the winding rate of a field line about the flux-tube axis, per unit 
length along the axis, is 

q = — = aJl ( kr ) r 17 \ 
LFF rB z krJa(kr) ' 

which varies in the radial direction but not in the vertical direction. At the 
flux-tube axis (r = 0), the winding rate is go,obs = lim r ^o 9 L ff = a /^- 

As the winding rate varies only in the radial direction, we can test the various 
fitting methods with error-free synthetic data generated with this magnetic field 
using the same approach used in Sections 4 and 5. To this end, in Equation (1) we 
set w h s ,i = fi<7 LFF (fi)- For synthetic data generated with this magnetic field it 
can be shown that the best-fit on-axis and average winding rates are not directly 
proportional to the respective true values. Therefore, we test the performance 
of the fitting methods over a range of winding rates, choosing a typical range of 
values as follows. For a flux loop of length I, the critical twist typically quoted for 
the onset of the kink instability is \ql\ w 2tt radians (e.g. Raadu, 1972; Hood and 
Priest, 1979; Einaudi and van Hoven, 1983; Mikic, Schnack, and van Hoven, 1990; 
Vclli, Hood, and Einaudi, 1990; Lionello et at, 1998; van der Linden and Hood, 
1998, 1999; Baty et at, 1998; Baty, 2001; Fan and Gibson, 2003, 2004; Torok, 
Kliem, and Titov, 2004). For the magnetic-flux tube described by Equation (13), 
we consider a vertical section of length SI = 4 Mm, taken to represent part of a 
larger magnetic-flux system. We construct twenty synthetic data sets with twenty 
equally spaced values for go,obs that yield partial twist values at the axis in the 
range — ir/2 < <?o,obs<^ < tt/2; the case where the winding rate at the axis is zero 
is not considered. We set m = 0.08 Mm -1 so that the magnetic- flux function 
decreases by approximately 27% from z = to z = SI = 4 Mm (at a fixed r). 
For this range of parameter values it can be confirmed that B z is single signed 
in the volume bounded by the flux surface * = an d < z < SI; we have 
restricted our attention to this volume to avoid locations where the winding rate 
for this field is infinite (corresponding to the zeroes of the Bessel function Jo, 
i.e. B z = 0). Synthetic measurements are generated at twenty equally-spaced 
observation points in the interval < < R, all of which are used for the 
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fitting, the average winding rate is calculated for the interval [0, R], and we set 
R = 2 Mm. For this choice of R the spacing between observation points (in the 
radial direction) is approximately 100 km, consistent with the spatial resolution 
provided by data from instruments such as the Solar Optical Telescope onboard 
Hinode (e.g. Kosugi et ai, 2007; Tsuneta et at, 2008). 




-0.4 -0.2 0.0 0.2 0.4 -0.4-0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 

9o.oi»( radians/Mm) g, ¥<ob ,( radians/Mm) qr . ob ,( radians/Mm) 



Figure 5. (a) The ratio of the inferred and true winding rates at the flux-tube axis [qg/qo,obs] 
as a function of the true winding rate at the axis [qo t, s ] for different fitting models. In 
these tests we use the linear force-free magnetic field (see Equation (17)) to generate twenty 
synthetic data sets with a range of winding rates go.obs with fixed m (see text). Each 
point represents the result from a single fitting experiment; the curves joining the points 
are only included as a guide. The purple curve is the best-fit case for a fitting model 
with q m (r) = q c , the blue curve is for q m (r) = q c (l — (r/R) 2 ) 2 , the green curve is for 
q m (r) = q c (l + 2(r/R) 2 - 3(r/R) 4 ), the gray curve is for q m (r) = Qi + Q2(r/R) 2 , and the 
orange curve is for q m {r) = Qi + Q2(r/R) 2 + Qs(r / R) 4 . (b) Same as (a) except for the 
corresponding ratio of the inferred and true average winding rates [q^ v / q av ,obs) as a function 
of the true average winding rate [q av b s ], with the average taken over the range used for the 
fitting, (c) The corresponding value of X 2 /{R 2 q^ obs ) (see Equation (1)) for the best-fit model 
as a function of q$ t, s . 

The results of this exercise (see Figure 5) are qualitatively similar to those 
discussed in Section 5. For example, we find that both the on-axis and average 
winding rates retrieved by the best-fit polynomial basis-function approximations 
are generally more accurate (and the \ 2 values are smaller) than those retrieved 
with fitting models with a fixed radial variation for the winding rate; although for 
the case shown in Figure 5 the magnitude of the discrepancy between the inferred 
and true winding rates is generally smaller than in Tables 1 and 2. We also find 
that the various fitting models which make fixed assumptions about the radial 
variation of the winding rate can produce different estimates for the winding 
rates given the same synthetic data, as in Table 2. Several additional points are 
worth noting: i) For the best-fit polynomial basis-function approximations the 
magnitude of the discrepancy tends to increase with increasing magnitude of the 
winding rate at a fixed value for m; the same is also true for the fitting model 
that assumes the winding rate is constant. This can be understood by referring 
to Figure 6(a) which shows how the winding-rate profiles (for a fixed value for 
m) vary with go.obs- For small values of go.obs the winding-rate profile is almost 
constant over the flux-tube interior, but as go, obs increases the change in the 
winding rate from the axis to the external boundary increases. This is because 
higher order terms in Taylor-series expansion for g LFF (r) about r = become 
more important as go, obs (or a) increases. Consequently, low-order polynomial 
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Figure 6. (a) Winding rate [<? LFF ] as a function of r (see Equation (17)) for a flux tube with 
m = 0.08 Mm -1 . The different curves correspond to four different values for the winding rate 
at the flux-tube axis: go.obs = a/2 = {0.05, 0.075, 0.1, 0.125}7r radians Mm -1 , (b) Winding rate 
[q LFF ] as a function of \I/ / S&o for a flux tube with m = 0.08 Mm -1 and a = 0.2-7T radians Mm -1 . 
The different curves correspond to five different heights: z = 0,1,2,3,4 Mm, with greater 
heights corresponding to curves with larger values for q LFF . 



basis-function approximations (and the fitting model with a uniform winding 
rate) are less accurate at large values of qo,obs- h) In other experiments with 
different values for m but the same range of values for go,obs we find that the 
magnitude of the discrepancy between the inferred and true winding rates tends 
to increase slightly as m increases. 

If the observed magnetic field varies in the z-direction, the results in Fig- 
ure 5 indicate that reasonable estimates for the winding rates can be obtained 
by fitting methods that infer the radial variation of the winding rate, for the 
plane where the magnetic field is measured, provided that the observed field is 
axisymmetric and has a flux-tube axis is straight and vertical. However, because 
the magnetic field varies in the z-direction the following considerations must be 
emphasised regarding the validity of the winding-rate estimates away from the 
plane where the magnetic field is measured. 

For an axisymmetric magnetic-flux tube, with an axis that is straight and 
vertical, and with a magnetic field that may vary in the z-direction, the winding 
rate at the flux-tube axis does not vary in the z-direction (see, e.g., Ferriz-Mas 
and Schussler, 1989, their Equation (4.15)). Thus, if an accurate estimate for 
the winding rate at the flux-tube axis can be retrieved at one height (by any 
method), it can be used to estimate the winding rate at the axis at other heights, 
provided that the observed flux tube is indeed axisymmetric with an axis that 
is straight and vertical. 

Away from the flux-tube axis, an axisymmetric magnetic-flux tube can gen- 
erally have a winding rate that varies in the z-direction (see, e.g., Browning and 
Priest, 1983, their Equation (5.7)). Therefore, estimates for the winding rate 
away from the flux-tube axis, and thus the average winding rate, retrieved with 
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fitting methods are generally valid only for the plane where the magnetic field is 
measured. The magnetic field described above is a special case where the winding 
rate has no z-dependence (see Equation (17)) and, therefore, the radial average 
of the winding rate computed over a fixed averaging interval is independent of z, 
but this is not expected in general for fields that vary in the vertical direction. 

From a different perspective, if the magnetic field varies in the vertical di- 
rection, the winding rate for a particular field line over a range of heights may 
be more important than the winding rate at a given radial location. For an 
axisymmetric magnetic-flux tube with a field that varies in the z-direction, the 
radial locations of the surfaces of constant ^ generally vary in the z-direction 
(away from the axis), see, e.g., Equation (13). Therefore, the winding rate for a 
particular field line at one height does not generally correspond to the winding 
rate of the same field line at some other height (away from the axis), see Fig- 
ure 6(b). Again, from this perspective if the observed magnetic field varies in the 
vertical direction, the winding-rate estimates retrieved by fitting methods away 
from the flux-tube axis are generally expected to be valid only for the plane 
where the magnetic field is measured. 

As mentioned above, the radial component of the magnetic field (see Equa- 
tion (14)) is non-zero (for m ^ 0). This suggests that the fitting model is not 
strictly appropriate for synthetic data constructed with this magnetic field. How- 
ever, the radial component is single- valued for a given value of r at fixed height, 
which would indicate to the observer that the magnetic field is axisymmetric (as 
assumed by the fitting model). The non-zero radial component of the magnetic 
field does not directly affect the winding-rate estimates. This is because the 
winding rate at a given position is defined as the angle that a field line rotates 
about the flux-tube axis per unit length along the axis and, thus, we are fitting for 
the ratio w = Bg/B z , which does not involve B r . If a non-zero radial component 
of the magnetic field is detected in the measurements this indicates that the 
magnetic field may be varying with height and, therefore, the caveats discussed 
above apply (provided that the observed flux tube is axisymmetric with an axis 
that is straight and vertical). 



B. Testing the Fitting Methods with Twisted Magnetic Fields that 
Vary in the Azimuthal Direction 

To show how the fitting methods perform when the observed flux tube has a 
magnetic field, and corresponding winding rate, that varies in the azimuthal 
direction we consider a magnetic field with: 

Be(r)=q c , aba rB cxp[-2(r/R) 2 ] , (18) 
B z (r, 9) = B cxp [-(r/R) 2 ] + B a \(r/R) 2 - (r/R) 3 ] sin (n a 6) , (19) 
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where q c obs is a winding rate, B and B a are magnetic-field strengths, R is the 
radius of the flux tube, and n a is an integer. At r = the magnetic field and cor- 
responding current density are both parallel to the z-direction, and the Lorentz 
force is zero. Away from r = 0, the azimuthal component of this magnetic 
field [Bg] varies only in the radial direction, whereas the vertical component 
[B z ] can vary in both the radial and azimuthal directions (for B a ^ and 
n a 7^ 0). This magnetic field is useful for testing the fitting methods because it 
has some features that are consistent with the fitting model {e.g. the flux-tube 
axis is straight and vertical, and the magnetic field does not vary in the vertical 
direction), which allows us to determine the influence of the variation of the 
magnetic field in both the azimuthal and radial directions on the winding-rate 
estimates retrieved with fitting methods. 

For B a = or n a = 0, the magnetic field described by Equations (18) and 
(19) is axisymmetric and the winding rate is q(r) = q c obs cxp(— (r/i?) 2 ), which 
has the same radial dependence as the synthetic data used in Table 2. On the 
other hand, for B a ^ and n a ^ 0, the winding rate of a magnetic- field line at 
the position (r, 9) about the flux-tube axis, per unit length along the axis, 



q = -L = q rabs exp [-(r/Rf] {l + (B a /B ) cxp [(r/i?) 2 ] (r / Rf [1 - (r/R)} sin (■ 
' o z 

(20) 

varies in both the r- and 0-directions. In either case, the winding rate at the 
flux-tube axis (r = 0) is q ,obs = 9 c , obs - 

To generate synthetic data, we sample the magnetic field described by Equa- 
tions (18) and (19) at a set of discrete locations on an x-y-plane, with a grid 
spacing of 100 km in both the x- and y-directions. We set B = 3.5 kG, 
R = 20 Mm, and n a — 25; these parameter values are chosen to be roughly 
consistent with a typical sunspot where azimuthal variations in the magnetic 
field are expected. We use all points with < R for the fitting, we calculate 
the radial average of the winding rate for the interval [0, R], and we choose the 
various parameter values for the synthetic data so that B z > within the flux 
tube. This magnetic field may not closely resemble the types of fields expected 
in sunspots, but we emphasise that the point of this exercise is to quantify how 
azimuthal variations in this field affect the estimates for the on-axis and average 
winding rates retrieved by fitting methods that use axisymmetric fitting models. 

Because the magnetic field and the winding rate vary in the azimuthal di- 
rection, the ratio u> bs,i derived from the synthetic measurements is not single- 
valued for a fixed value of r (see Figures 7 and 8). The parameter B a is related 
to the amplitude of the variation of the winding rate in the azimuthal direction 
(see Equation (20)), and we find that the scatter in the values of u» bs.i for a 
fixed value of r tends to increase as B a increases. These features indicate that 
the fitting model is not strictly appropriate for synthetic data constructed with 
this magnetic field. It may be possible to overcome this issue by extending the 
basis-function approximation approach of Section 5 to fit for both the azimuthal 
and radial dependence of the winding rate, but this is beyond the scope of 
the present investigation. Instead, we apply the fitting methods discussed so 
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Figure 7. The ratio, w = Bg/B z , as a function of r. The black points correspond to the 
ratio «J b s i derived from the synthetic measurements of a non-axisymmetric magnetic field 
(see Equations (18) and (19)) with q c obs = 0.1 radian Mm" 1 , B a /B = 1, R = 20 Mm, and 
n a = 25. The purple curve is the best-fit case for a fitting model with q m (r) = q c , the blue 
curve is for q m (r) = q c ((1 — (r/R) 2 ) 2 ), the green curve is for q m (r) = q c (l+2(r / R) 2 — 3(r/R) 4 ), 
the red curve is for q m (r) = q c exp(—(r/R) 2 ), the gray curve is for q m {r) = Qi + Q2(r / R) 2 , 
and the orange curve is for q m (r) = Qi +Q2(j'/R) 2 + Qi{f / R) 4 ', where R is the radius of the 
fitting model (here R = 20 Mm). 

far without modification (see, e.g., Figure 7). Because fitting methods that use 
axisymmetric fitting models cannot retrieve the azimuthal dependence of the 
winding rate, to test the efficacy of the fitting methods we use the azimuthally- 
averaged, radially-averaged winding rate of the observed magnetic field, along 
with the winding rate at the flux-tube axis. 

For synthetic data generated with this magnetic field, it can be shown that the 
best-fit on-axis and radially-averaged winding rates are directly proportional to 
the true on-axis and true azimuthally averaged, radially averaged winding rates, 
respectively. In Table 3 we show the ratio of the inferred and true winding rates at 
the flux-tube axis [<7o/<7o,obs] along with the ratio of the inferred radially averaged 
winding rate to the true azimuthally averaged, radially averaged winding rate 
[itv / lav ,obs] for three different values of B a /Bo. We include the axisymmetric 
case (with B a = 0) for reference and because the magnetic field for these tests 
is sampled differently to the case shown in Table 2. We find that the results of 
these tests are not strongly sensitive to the value of n a (although this depends on 
how the magnetic field is sampled). On the other hand, it is clear from Table 3 
that the results are sensitive to the value of B a /B . 

For the case with B a /B = 1, despite the presence of moderate azimuthal 
variations in the magnetic field and winding rate, we find the same general trend 
as in the cases that use axisymmetric synthetic data [B a — 0]. For example, we 
find that both the on-axis and average winding rates retrieved by the best-fit 
polynomial basis-function approximations are generally more accurate (and the 
X 2 values are smaller) than those retrieved with fitting models with a fixed 
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Table 3. Results for tests of several fitting models applied to synthetic data con- 
structed with the non-axisymmetric magnetic field described by Equations (18) and 
(19), with Bo = 3.5 kG, R = 20 Mm, and n a = 25, for three different values of 
B a /B . 



Fitting model «o/go,obs 9av/9av,obs X 2 /(-R 2 ^ obB ) 



B a /B = 





= 1c 


0.529 


0.708 


0.00885 


<? m ( r ) 


= q c (l-(r/R) 2 ) 2 


1.709 


1.220 


0.0511 


q m (r) 


= q c (l + 2(r/R) 2 -3(r/R) 4 ) 


0.562 


0.803 


0.0115 


1 m (r) 


= q c cM-(r/R) 2 ) 


1.0 


1.0 


0.0 


1 m (r) 


= Qi+Q2(r/R) 2 


0.902 


0.958 


0.000138 


Q m (r) 


= Qi+Q2(r/R) 2 + Q [i (r/R) 4 


0.988 


0.996 


9.586 xlO" 



Ba/Bo -- 


= 1 


Q m (r) = 




1 m ( r ) = 


q c (l-(r/R) 2 ) 2 


1 m ( r ) = 


q c (l + 2(r/R) 2 -3(r/R) 4 ) 


1 m ( r ) = 


q c e X p(-(r/R) 2 ) 


i m (r) = 


Qi+Q2(r/R) 2 


Q m (r) = 


Qi+Q2(r/R) 2 + Qa(r/R) 4 



0.539 


0.714 


0.0123 


1.741 


1.231 


0.0563 


0.575 


0.813 


0.0142 


1.019 


1.009 


0.00304 


0.923 


0.969 


0.00310 


0.983 


0.995 


0.00303 



B a /Bo = 2 





= Qa 


0.574 


0.736 


0.0274 


i m (r) 


= q c (l-(r/R) 2 ) 2 


1.852 


1.265 


0.0780 


i m ( r ) 


= q c (1 + 2(r/R) 2 - 3(r/R) 4 ) 


0.619 


0.845 


0.0265 


Q m (r) 


= q c c x p(-(r/R) 2 ) 


1.088 


1.040 


0.0166 


1 m (r) 


= Qi+Q2(r/R) 2 


0.997 


1.006 


0.0163 




= Qi+Q2(r/R) 2 + Q 3 (r/R) 4 


0.960 


0.990 


0.0163 
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radial variation for the winding rate. The exception is for the fitting model with 
q m (r) = q a exp(— (r/i?) 2 ), which retrieves accurate winding-rate estimates since 
it matches the winding-rate profile for the axisymmetric version of these data. 
We also find that the various fitting methods that make fixed assumptions about 
the radial variation of the winding rate can produce different estimates for the 
winding rates given the same synthetic data. 

For a given fitting model, the discrepancy between the inferred and true 
winding rates does not necessarily increase as the value of B a /B increases (as 
one may expect), although the magnitude of the % 2 values does generally increase 
(see Table 3). For the best-fit polynomial basis- function approximations applied 
to synthetic data with B a /B = 2, we find that the discrepancy between the 
inferred and true winding rates increases as nj, increases, contrary to the trend 
found at lower values of B a /Bo. Nevertheless, we generally find for the best- 
fit polynomial basis-function approximations that the discrepancy between the 
inferred and true winding rates is less than 10% for B a /Bo = 2 (see Table 3). 
Those discrepancies are quite small considering that the winding rate varies 
approximately by a factor of three from its minimum to maximum value for 
some values of r for B a /Bo = 2; for the case with B a /Bo = 1 the winding rate 
can vary by a factor of roughly two (see Figure 8(a)). We find similar results 
in other experiments with different functional forms for both the azimuthal and 
radial dependence of the magnetic field as those used in Equations (18) and (19) 
(results not shown). 




1 2 3 4 5 6 5 10 15 20 

0(radians) r(Mm) 



Figure 8. (a) Winding rate [q(r, 9)] as a function of 9 at fixed r for the non-axisymmetric 
magnetic field described in the text (see Equation (20)), with Bq = 3.5 kG, R = 20 Mm, 
Ba/Bo = 1, n a = 25, and q a obs = 0.1 radian Mm -1 . The red, green, blue curves correspond 
to r =5, 10, 15 Mm, respectively, (b) Winding rate [q(r, 9)] as a function of r at fixed 9 for the 
same flux tube as (a). The red, blue, green curves correspond to 9 =0 , 7r/50 , 37r/50 radians, 
respectively. 



C. Testing the Role of the Location of the Fitting-Model Axis 

Some methods that are commonly used to estimate the location of the axis of 
a magnetic-flux tube are: the peak of |J3|, the peak of \B Z \, the flux-weighted 
centroid, and the peak of the force-free parameter a z (e.g. Leamon et ai, 2003; 
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Lcka, Fan, and Barnes, 2005; Nandy et at, 2008). As discussed by Leka, Fan, 
and Barnes (2005) some of these methods do not consistently recover the correct 
location of the flux-tube axis; we have confirmed this finding using synthetic 
data constructed with a model for a toroidal magnetic- flux loop (e.g. Fan and 
Gibson, 2003, 2004, and see Section 6.1), but do not show the results for the 
sake of brevity. To demonstrate how the winding-rate estimates are affected by 
an error in the location of the fitting-model axis (in the absence of other sources 
of error), we generate synthetic data using a magnetic field with 

Be (r) = q c , oh jB cxp [-2(r/R) 2 } and B z (r) = B exp [-(r/R) 2 ] , (21) 

where q c oba is a winding rate, Bq is a magnetic-field strength, and R is the radius 
of the flux tube. To generate synthetic data we set q cobB = 0.1 radians Mm -1 , 
Bo = 3.5 kG, and R = 20 Mm, and sample the magnetic field at a set of discrete 
locations on an x-y-plane, with a grid spacing of 100 km in both the x- and 
y-dircctions. We introduce an error into the location of the flux-tube axis by 
computing synthetic data with the axis of the flux tube located at x = 4 Mm 
and y = 0. We apply the various fitting methods without modification assuming 
that the fitting-model axis is located at x = and y — 0. To estimate the winding 
rates we use only observation points that lie within 20 Mm of the fitting-model 
axis and have B z > 0. We calculate the average winding rate over the interval 
used for the fitting (i.e. < r < max(rj)). 

When there is an error in the location of the fitting-model axis we find that 
w is not single- valued for a fixed value of r (see Figure 9) , which indicates that 
the fitting model is not strictly appropriate for these synthetic data. Neverthe- 
less, it can be shown that the best-fit on-axis and average winding rates are 
directly proportional to the respective true values for this test. In Table 4 we 
provide the ratios gp/go.obs and <7av/9av,obs> along with the corresponding value 
for X 2 /(R 2 q 2 obs ) for the various best-fit models shown in Figure 9. The magnetic 
field described by Equations (21) is the same as for the axisymmetric case (with 
B a = 0) used in Appendix B; hence, for reference, the results in Table 3 for 
the case with B a — are those that would be retrieved if the fitting-model axis 
were correctly located for this synthetic data. Evidently, the accuracy of the 
winding-rate estimates retrieved by the fitting methods is affected by an error in 
the location of the fitting-model axis. However, for the better-performing cases 
(such as the polynomial basis-function approximations with rib — 3 and the case 
with the correct winding-rate profile) the magnitude of the relative discrepancy 
caused by the error in the location of the fitting-model axis is less than 10%. 
We find that the magnitude of the discrepancy in the inferred winding rates 
generally gets progressively larger as the error introduced into the axis location 
is increased. In experiments, using synthetic data generated with winding-rate 
profiles and magnetic fields different from Equation (21), we find that an error in 
the location of the fitting-model axis generally produces an error in the winding- 
rate estimates retrieved by the fitting methods of comparable magnitude to the 
case discussed here. 
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t — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — r 
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Figure 9. The ratio w = Bg /B z as a function of r (the radial distance from the fitting-model 
axis) for a test where an error is introduced into the location of the fitting-model axis, see Ap- 
pendix C for details. The black points correspond to the ratio w \> a i derived from the synthetic 
measurements. The purple curve is the best-fit case for a fitting model with q m (r) = q c , the blue 
curve is for q m (r) = q a ((1 — (r/R) 2 ) 2 ), the green curve is for q m (r) = q e (l+2(r / R) 2 — 3(r/R) 4 ), 
the red curve is for q m (r) = q c exp(— (r/R) 2 ), the gray curve is for <jf m (r) = Q± + Q2(r / R) 2 , 
and the orange curve is for q m (r) = Qi + Qziyl R) 2 + Q:j(r / R) 4 , where R is the radius of the 
fitting model (here R = 20 Mm). 

Table 4. Summary of results for a test where an error is introduced into the location 
of the fitting-model axis, see Appendix C for details. 

Fitting model ^/^O.obs gav/9av,obs X 2 /{R 2 <ll _ obs ) 



9™ 


(r) 


= la 




0.539 


0.722 


0.00841 


1m 


(r) 


= q c {l-(r/R) 2 ) 2 




1.573 


1.124 


0.0445 


Irr, 


(r) 


= q c (l + 2(r/R) 2 - 


3(r/i?) 4 ) 


0.524 


0.748 


0.0121 


1-m 


(r) 


= q c exp(-(r/R) 2 ) 




0.957 


0.957 


0.00323 


1m 


(r) 


= Qi+Q 2 (r/R) 2 




0.820 


0.893 


0.00323 


% 


(r) 


= Qi+Q2(r/R) 2 4 


-Q S (r/R,) 4 


0.926 


0.940 


0.00301 
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